Negative specific heat in a Lennard- Jones-like 
gas with long-range interactions 

Ernesto P. Borges a ' 5 * Constantino Tsallis" 

a Centro Brasileiro de Pesquisas Fisicas, 
R. Dr. Xavier Sigaud 150, 22290-180 Rio de Janeiro-RJ, Brazil 
b Escola Politecnica, Universidade Federal da Bahia, 
R. Aristides Novis, 2, 40210-630 Salvador-BA, Brazil 



Abstract 

We study, through molecular dynamics, a conservative two-dimen- 
sional Lennard- Jones-like gas (with attractive potential oc r~ a ). We 
consider the effect of the range index a of interactions, number of 
particles, total energy and particle density. We detect negative specific 
heat when the interactions become long-ranged (0 < a/d < 1). 
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Non-Gaussian equilibrium in long-range Hamiltonian systems has recently 
been reported 0. The authors studied a model of classical planar spins with 
infinite range interaction, and showed that before going to the Boltzmann- 
Gibbs equilibrium state, the system reaches a metastable state with a non- 
Maxwellian distribution of velocities. The duration of this metastable state 
appears to diverge with the number of particles, i.e., in the thermodynamical 
limit, the system indefinitely stays trapped in a non-Boltzmann-Gibbs state. 
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In this paper, we investigate similar anomalies with a different long- 
range Hamiltonian model. We consider a two-dimensional gas confined in 
a box of linear length L, defined by the Hamiltonian Ti = K + V/ N + 
V wa ii s , where K = S^ 1 p|/(2m), the interacting potential is given by V = 



Ef<^(r !3 ) with @ v( nj ) = C a [(cr/r tJ ) 12 -(a/r tJ ) a \ (0 < a < 12) and 

C a = e (12 12 /a a ) 1 '( 12 ~ a > /(12 — a), is the distance between the centers of 
particles i and j, a characterizes the diameter of a particle, m is its mass and e 
is the energy scale. We adopt appropriate units in which a — e = m = 1. C a 
is a factor that makes the depth of the potential well the same Va. V wa ii s is the 
repulsive potential of the walls, which assures that the particles are confined 
within the box. We adopt soft walls, namely V wa a s = Y,ti=i J2f=i(^/ r iw) 12 , 
where w is the index for the walls, and r iw is the distance between particle i 
and wall w. The nonextensive scaling parameter N = 1 + d dr r d_1_Qf 
(see ||) is convenient to make the Hamiltonian formally extensive \/a/d (d 
is the spatial dimensionality of the system). The model is a variation of the 
one introduced in J4|, and the particular case a = 6 recovers the celebrated 
Lennard- Jones model for dilute nonpolar gases. The variable involved in 
the range of interactions (r^) is also the dynamical variable, while in the clas- 
sical spins model, the range of interactions is controlled by the (fixed) position 
of the particles in the lattice, and not by the angles, which are the dynam- 
ical variables. We integrated the equations of motion using the forth order 
symplectic algorithm || with time steps At = 0.005-0.02, which yielded 
a maximum relative error in the total energy conservation AE/E ~ 10 -5 . 
We adopted as initial conditions for positions (a) all particles at a triangu- 
lar lattice (note that this geometry corresponds to the global minimum for 
a two-dimensional L 2 system), and (b) particles randomly distributed, and 
a water bag distribution for initial velocities. We considered a long-range 
system with a/d = 0.5, and a short-range a/d — 3, for comparison. 

Fig. |l] shows caloric curves for short- and long-range systems with 100 
particles. Each point of the caloric curve is obtained by a time average after 
a transient of the order of 10 2 -10 3 time units (see Inset), and by an average 
of a few (up to 10) events with different initial conditions. A negative specific 
heat region is exhibited for a/d = 0.5, while this effect does not appear in 
the short-range case. 

Data collapse in the abscissa of the caloric curve (Fig. 0) for the long-range 
case is achieved when the total energy scales as E/N 1+x , with x = 0.0876. 
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Such small value of x is not easily distinguishable (numerically speaking) from 
a possible logarithmic correction hidden in the power-law scaling introduced 
by N. Such effect was not detected in the classical spins model []TJ, and 
could be due to the boundary conditions. Data collapse in the ordinate of 
the caloric curve may be roughly achieved by scaling the maxima of the curves 
(for energies below the critical) according to T max (N) ~ T max (oo) — a/N y 
(T = (K)/N) with y = 0.038, a value once again very small. 

Fig. [| displays the effect of the particle density in the caloric curve. The 
negative specific heat region increases with decreasing density, when long- 
range forces are present. Highly dense systems don't exhibit such anomaly. 

Finally, Fig. ^ shows the probability distribution function (pdf ) of veloci- 
ties for the long-range (a/d = 0.5) system, with N/L 2 = 10~ 3 and E/N = 0.7 
(a state in the negative specific heat region). As N increases, the curves seem 
to departure from a Maxwellian, exhibiting a slightly positive curvature. A 
non-Maxwellian pdf was recently found in the classical spins model PJ, that 
might be understood within nonextensive statistical mechanics || [7], H, !]• 
Such framework may also be the proper one to handle classical gases with 
long-range interactions, but only further results would enable a clear conclu- 
sion. 
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support and hospitality. We also acknowledge CAPES, FAPERJ, CNPq, 
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Figure 1: Caloric curves for N = 100 and N/L 2 = 10 2 : a/d = 3 (squares), 
a/d = 0.5 (circles). Solid lines are guides to the eyes. Inset: Typical time evolution 
of (K)/N. 



5 



0.4 



0.3 

<K>/N 

0.2 



0.1 



0.0 



1.0 

ca 

&0.9 

o> 
o 

Z 0.8 



Q. 

03 
3 
O 



lu°0.7 



o/d = 0.5 a, 
N/L 2 = 10" 2 . 



10 100 N 1000 ^ 



N= 50 • 

100 ■ 

200 ♦ 

400 a 







0.2 



0.4 



E/N 



1+x 



0.6 



0.8 



Figure 2: Effect of the number of particles in the caloric curve, for the long-range 
system (a/d = 0.5, N/L 2 = 1(T 2 ). N = 50, 100, 200, 400. Inset: the energy per 
particle at the cusp point as a function of N (log-log scale), slope x = 0.0876. 
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Figure 4: Probability distribution function for the long-range case (a/d = 0. 
E/N = 0.7 (below the critical energy per particle) and N/L 2 = 1(T 3 . N = 81 
1600, 2200. We have taken, in the abscissa, the values of both p x and p y . 
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